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Abstract. Based on a simple, but fairly successful, model 
of the chemical evolution of the Milky Way disk, we study 
the evolution of the abundances of the elements He, C, N, 
O, Ne, Mg, Al, Si, S, Ar and Fe. We use metallicity depen- 
dent yields for massive stars with and without mass loss. 
We find that most observed abundance profiles are cor- 
rectly reproduced by massive star yields, but C and N re- 
quire supplementary sources. We argue that massive, mass 
losing stars can totally account for the abundance profile 
of C, while intermediate mass stars are the main source 
of N; in both cases, some conflict with corresponding data 
on extragalactic HII regions arises, at least if current ob- 
servations in the Galaxy are taken at face value. The ob- 
served behaviour of Al is marginally compatible with cur- 
rent massive star yields, which probably overestimate the 
"odd-even" effect. We also find that the adopted "inside- 
out" formation scheme for the Milky Way disk produce 
abundance profiles steeper in the past. The corresponding 
abundance scatter is smaller in the inner disk than in the 
outer regions for a given interval of Galactic age. 
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1. Introduction 

Abundance gradients constitute one of the most important 
observational constraints for models of the evolution of the 
Milky Way disk. The existence of such gradients is now 
well established, through radio and optical observations of 
HII regions, stars and planetary nebulae (see Henry and 
Worthcy 1999 for a review of the abundance profiles in 
the Milky Way as well as in external galaxies). An average 
gradient of dlog(X/H)/dR~ —0.06 dex kpc -1 is observed 
in the Milky Way for O, S, Ne and Ar, while the gradients 
for C and N are, perhaps, slightly larger (Smartt 2000). 
Observations of open clusters support the above picture 
(Friel 1995, 1999; Carraro et al. 1998), albeit with large 
uncertainties. 



The existence of those gradients offers the opportunity 
to test theories of disk evolution and stellar nucleosynthe- 
sis. Indeed, the magnitude of the observed gradients in the 
Milky Way disk is rather large: metal abundances in the 
inner disk are larger than those in the outer disk by a fac- 
tor of ~10. This suggests that the role of the Galactic bar 
in inducing large scale radial mixing and therefore flaten- 
ning the abundance profile (e.g. Friedli and Benz 1995) 
has been rather limited; alternatively, the Galactic bar is 
too young (<1 Gyr) to have brought any important mod- 
ifications to the gaseous and abundance profiles of the 
disk. In any case, the magnitude of the observed gradi- 
ents puts important constraints on disk models invoking 
radial inflows (see e.g. Portinari and Chiosi 2000 and ref- 
erences therein). Even in the case of simpler models (with 
no radial inflows), the abundance gradients constrain the 
various parameters, like the local timescales of star forma- 
tion and infall (e.g. Matteucci and Francois 1989, Prantzos 
and Aubert 1995, Molla et al. 1997) or any variation of 
the stellar Initial Mass Function properties with metallic- 
ity (e.g. Chiappini et al. 2000). On the other hand, profiles 
of abundance ratios across the disk constrain the nature of 
various elements and isotopes (i.e. "secondaries" vs. "pri- 
maries", "odd-Z" vs. "even-Z") and their nucleosynthesis 
sites (e.g. Prantzos et al. 1996). 

A large number of chemical evolution models has been 
developped in the 90ies, aiming to explore one or more of 
the above issues, either with radial inflows (e.g. Gotz and 
Koppen 1992, Chamcham and Tayler 1994, Tsujimoto et 
al. 1995, Firmani et al. 1996, Thon and Meusinger 1998, 
Portinari and Chiosi 2000) or without radial flows (Fer- 
rini et al. 1994, Prantzos and Aubert 1995, Giovagnoli 
and Tosi 1995, Carigi 1996, Chiapinni et al. 1997, Boissier 
and Prantzos 1999, Chang et al. 1999). Some promising 
chemodynamical evolutionary models for the Milky Way 
have also been developed by a few groups (Steinmetz and 
Miiller 1994, Samland et al. 1997, Berczik 1999). For the 
simplest of those models (namely those with no radial 
flows), some convergence has been recently reached be- 
tween the various groups, at least concerning the basic 
ingredients (e.g. Tosi 2000): i) necessity of substantial in- 



2 



J.L. Hou, N. Prantzos and S. Boissier: Abundance gradients and their evolution in the Milky Way disk 



fall and of radially varying timescales for the infall and the 
star formation, and ii) no need for varying IMF or strong 
galactic winds. 

Despite this agreement, some important differences 
still exist between the various models. The most impor- 
tant concerns the history of the abundance profiles: were 
they steeper or flatter in the past? The former is suggested 
by models of e.g. Prantzos and Aubert (1995), Molla et al. 
(1997), Allen et al. (1998), Boissier and Prantzos (1999), 
while the latter is supported by models of Tosi (1988) and 
Chiappini et al. (1997). The situation is not settled ob- 
servationally either. Estimated ages of open clusters and 
planetary nebulae of various types span a large fraction of 
the age of the Galaxy. Observations of the abundances of 
those objects across the Milky Way disk could, in princi- 
ple, provide some information on the past history of the 
abundance gradients. In practice, however, observational 
uncertainties are too large at present to allow conclusions 
on that matter (e.g. Carraro et al. 1998, Maciel 1999). 

Several observational studies on the abundance pro- 
files of the Milky Way disk were performed in the late 
90ies. They concerned a variety of elements (He, C, N, 
O, Ne, Mg, Al, Si, S, Ar and Fe), observed in HII re- 
gions (Simpson et al. 1995, Rudolph et al. 1997, Afflerbach 
et al. 1997, Deharveng et al. 2000), B-stars (Smartt and 
Rolleston 1997, Gummersbach et al. 1998) and planetary 
nebulae (Maciel and Koppen 1994, Maciel and Quireza 
1999). Although they leave the question of the history of 
abundance profiles unsettled, these studies allow for some 
important tests of stellar nucleosynthesis theories, since 
they provide absolute abundances, and abundance ratios 
along the Milky Way disk. 

In a previous work (Boissier and Prantzos 1999, here- 
fater BP99) we developed a simple model for the chemical 
and spectrophotometric evolution of the Milky Way disk. 
With few assumptions (basically concerning the radial 
variation of the infall and star formation timescales), the 
model reproduces all the major obsevational constraints 
for the Galaxy: total amounts of gas and stars, super- 
novae rates, radial profiles of gas, stars, star formation 
rate and oxygen abundance, luminosities and scalelengths 
in various wavelength bands. The success of the model 
does not guarantee, of course, its correctness or its unique- 
ness. However, it offers a sound basis for a more detailed 
exploration of issues related to the chemical evolution of 
the Milky Way. Two such issues are explored here: 

i) how successful are currently available stellar yields 
in reproducing the wealth of recently observed abundance 
profiles in the Galaxy? are Wolf-Rayet (massive, mass los- 
ing) stars or intermediate mass stars (IMS) important in 
shaping the abundance profiles, and for which elements? 

ii) is the evolution of abundance gradients predicted by 
the model compatible with observations of available trac- 
ers? are there any new, potentially testable, predictions of 
the model in that respect? 



The plan of the paper is as follows: In Sect. 2, we 
present briefly the basic ingredients and the underlying 
assumptions of the BP99 model; in particulrar, we present 
in some detail the metallicity dependent yields of Woosley 
and Weaver (1995) and Maeder (1992), on which much of 
our work is based. Model results for abundance gradients 
of He, C, N, O, Ne, Mg, Al, Si, S, Ar and Fe and their 
evolution are given in Sect. 3. Sect. 4 constitutes the main 
body of the paper. In sect. 4.1 we present a compilation 
of observational data from HII regions, B-stars and plane- 
tary nebulae in the Milky Way disk. Detailed comparison 
between model predictions and current observed disk pro- 
files is made in Sect. 4.2. In Sect. 4.3, we further compare 
the behaviour of model profiles of abundance ratios to 
observations and draw some conclusions on the adopted 
yields. Then, based on the abundance data for objects of 
different ages we discuss the evolution of abundance gradi- 
ents for elements O, Ne, S and Ar (Sect. 4.4). The results 
are summarized in Sect. 5. 

2. Model 

The adopted model for the chemical evolution of the Milky 
Way disk is described in detail in BP99. In Sect. 2.1 we 
briefly recall the main features of the model. In Sect. 2.2 
we present in some more detail the only novel ingredient 
with respect to BP99, namely the metallicity dependent 
yields of Woosley and Weaver (1995, hereafter WW95) for 
intermediate mass elements and of Maeder (1992, here- 
after M92) for He, C, N, O (in BP99 only yields for stars 
of solar metallicity are used). 

2.1. Description of the model 

The galactic disk is considered as an ensemble of con- 
centric, independently evolving rings, progressively built 
up by infall of primordial composition. The assumption 
of infall is traditionally based upon the need to explain 
the locally observed metallicity distribution of long-lived 
stars, which cannot be explained by the simple "closed- 
box" model (leading to the well-known "G-dwarf prob- 
lem"). However, the recent work of Blitz et al. (1999) gives 
observational support to this idea, showing that the Milky 
Way and M31 are currently accreting substantial amounts 
of gas (~ 1 M /yr) in the form of high velocity clouds of 
low metallicity. 

The infall rate is assumed to be exponentially decreas- 
ing in time, i.e. 

f(t,R) = A(R) e -'M*) (1) 

with a characteristic timescale t(Rq) = 7 Gyr in the so- 
lar neighborhood (Rq = 8kpc), in order to reproduce the 
local G-dwarf metallicity distribution. t(R) is assumed to 
increase outwards, from r(R = 2 kpc) = 1 Gyr to t(R = 17 
kpc) = 12 Gyr. This radial dependence of the timescale 
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of the infall rate f(R) is simulating the inside-out for- 
mation of galactic disks and, combined with the adopted 
SFR ^(i?) (Eq. 3) allows to reproduce the observed cur- 
rent profiles of gas, oxygen abundance and SFR in the 
Milky Way disk (see BP99 and Sects. 2.3 and 4 below). 
The coefficient A(R) is obtained by the requirement that 
at time T = 13.5 Gyr the current mass profile of the disk 
is obtained, i.e. 

f f(t,R) = E(fl) (2) 
Jo 

with oc e - R / R o anc i a scalelength Rq = 2.6 kpc for 

the Milky Way disk. 

The chemical evolution of each zone is followed by solv- 
ing the appropriate set of integro-diffcrential equations, 
without the Instantaneous Recycling Approximation. The 
adopted stellar Initial Mass Function (IMF) is a multi- 
slope power-law between 0.1 M Q and 100 M© from the 
work of Kroupa ct al. (1993), leading to a Return Frac- 
tion R = 0.32. 

The star formation rate (SFR) is locally given by a 
Schmidt-type law, i.e. proportional to some power of the 
gas surface density T, g : ^ oc S^ 5 , according to the obser- 
vations of Kennicutt (1998). It varies with galactocentric 
radius R, as: 

$(t,R) = a E^t.iZ) 1 - 5 V(R) R- 1 (3) 

where V(R) is the circular velocity at radius R. This radial 
dependence of the SFR is suggested by the theory of star 
formation induced by density waves in spiral galaxies (e.g. 
Wyse and Silk 1989). Since V(R) ~ constant in the largest 
part of the disk, this is equivalent to ^(R) oc T, g (R) 15 R~ 1 . 
The efficiency a of the SFR in Eq. 3 is fixed by the re- 
quirement that the local gas fraction <j g (R ) ~ 0.2, is re- 
produced at T = 13.5 Gyr. 

We assume that the "rings" of the disk are evolving in- 
dependently from one another. This (ovcr)simplification 
ignores in general the possibility of radial inflows in 
gaseous disks, resulting e.g. by viscosity or by infalling 
gas with specific angular momentum different from the 
one of the underlying disk; in both cases, the resulting 
redistribution of angular momentum leads to radial mass 
flows. The magnitude of the effect is difficult to evaluate, 
because of our poor understanding of viscosity and our ig- 
norance of the kinematics of the infalling gas. Models with 
radial inflows have been explored in the past (Mayor and 
Vigroux 1981; Lacey and Fall 1985; Clarke 1989; Cham- 
cham and Tayler 1994). It turns out that for some combi- 
nations of the parameters of infall, radial inflow and SFR, 
acceptable solutions are obtained, i.e. the current radial 
profiles of various quantities are successfully reproduced 
(see, e.g. Portinari and Chiosi 2000 for a recent overview of 
the problem). However, at the present stage of our knowl- 
edge introduction of radial inflows in the models would 



imply more free parameters than observables. For sim- 
plicity reasons we stick to the model of "independently 
evolving rings" for the Milky Way disk. 

2.2. Yields of massive stars 

An important ingredient in our study of abundance gradi- 
ents are the stellar yields of various elements. Most of the 
intermediate mass elements studied here are produced by 
massive stars, with the exception of some CNO isotopes 
that are also produced by intermediate mass stars. We 
consider no yields from intermediate mass stars in this 
work; in the line of Goswami and Prantzos (2000), con- 
cerning the evolution of the halo+local disk, our explicit 
purpose is to check to what extent massive stars can ac- 
count for observations of intermediate mass elements and 
for which elements the contribution of intermediate mass 
stars is mandatory. 

We use the metallicity dependant yields of WW95, 
which are given for stars of mass M = 12, 13, 15, 18, 
20, 22, 25, 30, 35 and 40 M and metallicities Z/Z = 0, 
10~ 4 , 10~ 2 , lO" 1 and 1. In Fig. 1 we present the WW95 
yields, folded with the Kroupa et al. (1993) IMF. They are 
presented as overproduction factors, i.e. the yields (ejected 
mass of a given element) are divided by the mass of that 
element initially present in the part of the star that is 
finally ejected: 

lMi^(M)^M)dM 
X e>i (M - M R ) *(M) dM 

where: $(M) is the IMF, Ml and M2 the lower and up- 
per mass limits of the stellar models (12 M Q and 40 M Q , 
respectively). Yi(M) are the individual stellar yields and 
Mr the mass of the stellar remnant. Adopting Xq^ in Eq. 
(4) creates a slight inconsistency with the definition of the 
overpoduction factor given above, but it allows to visual- 
ize the effects of metallicity in the yields of secondary and 
odd-Z elements. 

As can be seen from Fig. 1: i) most of the intermediate 
mass elements are nicely co-produced (within a factor of 
2) by solar metallicity stars; ii) the "odd-even effect" , fa- 
voring the production of odd- nuclei at high metallicities, 
is clearly present; iii) the yields of Ne and Mg show, curi- 
ously, some dependence on metallicity (not as large as the 
one of the odd-elements Na and Al, but still enough to lead 
to some interesting abundance patterns, as we shall see in 
Sect. 3); iv) He, C, N, Sc, V and Ti are underproduced 
relative to Oxygen. He, C and N clearly require another 
source (intermediate mass stars and/or Wolf-Rayet stars, 
see Prantzos et al. 1994 and Sect. 4.2), while the situation 
is less clear for the elements, Sc, V and Ti (see Goswami 
and Prantzos 2000). 

The calculations of WW95 did not consider any mass 
loss during stellar evolution. Thus, they probably underes- 
timated the yields of several elements that are expelled by 
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Fig. 1. Average overproduction factors (over a Kroupa 
et al. (1993) IMF, see Eq. 4) of the yields of Woosley 
and Weaver (1995) for 3 different initial stellar metallici- 
ties, covering reasonably well the metallicity evolution of 
the Milky Way disk. The solid horizontal line is placed 
at F oxygen and the two dotted horizontal lines at half and 
twice that value, respectively. Most of the intermediate 
mass elements are nicely co-produced (within a factor of 
two). The "odd-even effect" is clearly seen (e.g. in the 
cases of Na, Al, P, etc.); notice the small (but significant) 
metallicity dependence of the Ne and Mg yields. N behaves 
as a pure "secondary" element. The elements He, C, N, Li 
and Be obviously require another production site. 

the intense winds of massive stars, i.e. He, N and C. The 
effect of stellar winds is stronger when the stellar metallic- 
ity is larger. Maeder (1992) found that stars with M >30 
M Q and Z>0.1 Z© eject considerably larger amounts of 
He, N and C during their lifetime than their lower metal- 
licity counterparts; for that reason, less matter is left in the 
He-core to be processed into Oxygen. At lower metallici- 
ties, the effect of stellar winds is negligible and stars of all 
masses evolve almost at constant mass. Elements heavier 
than Oxygen are produced in the subsequent, very rapid, 
stages of stellar evolution (after core He exhaustion) and 
their yields are not directly affected by the intensity of the 
mass loss. However, the structure of the stellar core may 
be affected by the loss of mass and this may also affect the 
final yields of heavy elements (e.g. Woosley, Langer and 
Weaver 1993). 

In Fig. 2 we present the yields of M92 for He, N, 
C and O as a function of stellar mass; they are given 
for two metallicities (Z/Z Q = 0.05 and 1, respectively) 
and are compared to the corresponding yields of WW95. 
The aforementionned effect of metallicity-dependent stel- 
lar winds on the yields of stars with M>30 M Q is clearly 
seen. In Sect. 3.4 we shall explore the effect of those yields 
on the abundance gradients in the disk. 



i WW95 □ M92 
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Fig. 2. Massive star yields of He, C, N and O for differ- 
ent initial metallicities, according to WW95 (filled trian- 
gles) and M92 (open squares); Solid curves: Z = Zq, dotted 
curves: Z = 0.1 Z Q (for WW95), dashed curves: Z = 0.05 
Zq (for M92). There is good general agreement between 
the two calculations for stars with initial metallicity Z = 
Zq. He, C and O are primary elements, i.e. their yields 
are independant of the initial stellar metallicity; the dif- 
ferences between the two sets of M92 yields for stars with 
M>30 M are due to the metallicity dependent stellar 
winds. Nitrogen is produced as a secondary element in 
both WW95 and M92. 



To account for the additional source of Fe-peak ele- 
ments, required to explain the observed decline of O/Fe 
abundance ratio in the disk (e.g. Goswami and Prantzos 
2000) , we utilise the recent yields of SNIa from the explod- 
ing Chandrashekhar-mass CO white dwarf models W7 and 
W70 of Iwamoto et al. (1999). These are updated versions 
of the original W7 model of Thielemann et al. (1986), cal- 
culated for metallicities Z = Z Q (W7) and Z = (W70), 
respectively. In this model, the deflagration is starting in 
the centre of an accreting white dwarf, burns ~ half of 
the stellar material in Nuclear Statistical Equilibrium and 
produces ~ 0.7 M of 56 Fe ( in the form of 56 Ni). These 
SNIa models lead to an oveproduction of Ni, but the evo- 
lution of this element will no be considered here. 

It should be emphasised that the evolution of the SNIa 
rate is not well known, and hardly constrained by obser- 
vations. For the purpose of this work, we shall adopt the 
formalism of Matteucci and Greggio (1986) for the rate of 
SNIa, adjusting it as to have them appearing locally after 
the first Gyr, i.e. at a time when [Fe/H]~ —1 in the solar 



J.L. Hou, N. Prantzos and S. Boissier: Abundance gradients and their evolution in the Milky Way disk 



5 




4 8 12 16 



Radius (kpc) 

Fig. 3. Results of the chemical evolution model for the 
Milky Way disk and comparison to observations. In all 
panels, results are shown at three different epochs (dot- 
ted curves: 1 Gyr, dashed curves: 5 Gyr, solid curves: 
13.5 Gyr). The latter are to be compared to observations 
of present-day profiles in the Milky Way disk, shown as 
shaded regions in first three panels and within error bars 
in the lower panel. In all panels, the error bar at Rq = 
8 kpc indicates observed quantities in the solar neighbor- 
hood. References for data are given in BP99. 

neighborhood. However, use of that same formalism along 
the disk will lead to different O/Fe abundance ratios at 
T = 13.5 Gyr as we shall see in Sect. 3, i.e. the final Fe 
gradient will be different from the one of oxygen (see also 
Prantzos and Aubert 1995). 

2. 3. Results for the Milky Way disk 

As described in detail in BP99, the simple model pre- 
sented in Sect. 2.1 can readily account for the evolution of 
the solar neighborhood, reproducing quite successfully the 
main observational constraints (age-metallicity relation- 
ship, metallicity distribution of long-lived F-stars, current 



local surface densities of stars, gas, star formation and su- 
pernova rates). Also, in Goswami and Prantzos (2000) it is 
shown that the use of the WW95 yields for massive stars 
and of the Iwamoto et al. (1999) yields for SNIa leads to a 
successful agreement between the gaseous composition of 
the model at an age of 9 Gyr and the observed solar one. A 
few exceptions concern the elements He, C and N (which 
are underproduced) and Ni (which is overproduced, be- 
cause of the adopted SNIa yields). 

The model predictions for the disk are equally success- 
ful, at least to a first order. Indeed, the adopted combina- 
tion of SFR (Eq. 3) and infall rate (Eq. 1) leads to final 
profiles of gas and SFR that are in fair agreeement with 
the observed ones, as can be seen in Fig. 3. The stellar 
profile is also in agreement with observations, but it is 
essentially determined by the boundary condition of Eq. 

2. However, the adopted inside-out formation scheme of 
the disk leads naturally to different scalelengths in the B- 
band (reflecting mostly the SFR profile in the past ~ 1 
Gyr) and the K-band (reflecting the total stellar popula- 
tion, cumulated over T = 13.5 Gyr). As shown in BP99, 
the corresponding scale-lengths (~ 4 kpc in the B-band 
and ~ 2.6 kpc in the K-band, respectively) are in fair 
agreement with observations. Moreover, the model also 
reproduces reasonably well the total current SFR and su- 
pernova rates as well as the total luminosities in various 
wavelenght bands. This is a rather encouraging success, 
since the number of the constraints is much larger than 
the number of the parameters. The success of this simple 
model encourages us to use it for a thorough study of the 
various abundance gradients in the Milky Way disk. 

3. Model results for abundance gradients 

We calculated the evolution of the abundances of all el- 
ements between H and Zn with the WW95 metallicity 
dependent yields in all the zones of our model disk. In 
Fig. 4 we present the results concerning all elements with 
measured abundance gradients in the Milky Way. We no- 
tice that, since we did not include yields from IMS in our 
calculation, our results for He, C and N represent rather 
lower limits. We shall discuss the comparison to observa- 
tions in the next section, where we shall also explore the 
role of massive, mass losing stars to the abundance pro- 
files of those elements. Here we focus on the model results, 
which may be summarised as follows: 

i) Final values (at T = 13.5 Gyr) of the abundances 
at Ro = 8 kpc are ~ solar for O, Al, Si, S, Ar and Fe; 
they are slightly lower than solar for Ne and Mg; and 
they are considerably below solar for C and N. It may 
appear surprising that Ne and Mg do not reach their solar 
values, but this is a consequence of the adopted metallicity 
dependent yields of WW95 (Fig. 1): the average (i.e. over 
the disk's age) overproduction factors of Ne and Mg are 
lower than the one of e.g. O. It is not clear whether such 
a dependence is physical or just an artifact of the WW95 
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Fig. 4. Evolution of abundance profiles in the Milky Way 
disk, in the framework of the model presented in Sect. 2 
and with the WW95 metallicity dependent yields. Curves 
show results at 2, 4, 9 and 13.5 Gyr, respectively (as in- 
dicated in the lower right panel). Because of the metallic- 
ity dependence of the yields, some elements (N, Al) have 
steeper abundance profiles than e.g. O, at all times. This is 
also the case for Fe, but the reason is its production mostly 
by SNIa. C and N are underproduced with respect to their 
solar system values, because yields from intermediate mass 
stars are neglected in this calculation. Metallicity depen- 
dent massive star yields from M92 considerably improve 
the situation concerning carbon (see Fig. 7), but not ni- 
trogen. Solar abundances (0) were taken from Grevesse 
& Sauval (1998). 

models. As for C and N, it is clear that massive stars with 
no mass loss cannot be the main sources of those elements 
(see also Sect. 4). 

ii) The most prominent feature of the model is the pre- 
diction that abundance gradients flatten with time. This 
is a generic feature of all models forming the galactic disk 
"inside-out" . Indeed, in that case, there is a rapid increase 
of the metal abundance at early times in the inner disk, 
leading to a steep abundance gradient. As time goes on, 
star formation "migrates" to the outer disk, producing 
metals there and flattening the abundance gradient. 

hi) The final abundance profile (T = 13.5 Gyr) is, in 
general, flatter in the inner disk. As already described in 
Prantzos and Aubert (1995) this is due to the fact that 
in those regions the large populations of low-mass, long- 
lived stars that are formed early on in galactic history 
reject a lot of metal-poor gas at the end of their evolu- 



Fig. 5. Evolution of abundance gradients ( expressed in 
dex/kpc ) in the Milky Way disk (range: 4—14 kpc from 
the centre), in the framework of the model presented in 
Sect. 2 and with the WW95 metallicity dependent yields. 
Because of the inside-out formation scheme adopted here, 
gradients were steeper in the past. Because of the metal- 
licity dependence of the yields of WW95, some elements, 
like N or Mg, have larger abundance gradients (in absolute 
value) than O. This is also the case for Fe, but the reason 
is its production mostly by SNIa. 



tion, which dilutes the metal abundances; this effect is 
absent in the outer regions, where there are not very old 
stellar populations. Notice that the flattening seems to be 
more important in the case of N, Ne, Mg and Al. These 
elements show some metallicity dependence in their yields 
(at least according to WW95) which is difficult to under- 
stand in the case of Ne and Mg, but expected in the case 
of the secondary N and of the odd-Z Al. Since WW95 
give yields only up to stellar metallicities of Z = Z , for 
higher metallicities we use their Z = Z Q yields. This means 
that in the inner disk, where metallicities higher than so- 
lar are reached, we underestimate the production of any 
element with metallicity dependent yield. If appropriate 
yields were used, the Al and N profiles in the inner disk 
would be steeper, not flatter, than the one of O. 

The magnitude of the current abundance gradients in 
the Milky Way disk is one of the most important con- 
straints in the models of the evolution of our Galaxy. Most 
of the proposed models reproduce it fairly well (e.g. Tosi 
2000 and references therein), at least when no radial in- 
flows are included. However, equally important is the ques- 
tion of the evolution of those gradients and, in particular, 
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whether they flatten or steepen with time. We shall con- 
front our models to the data and to other theoretical works 
in the next section. Here, we present the evolution of our 
model gradients in the 4—14 kpc region for a few selected 
elements (Fig. 5). All gradients were systematically larger 
in the past. The gradient of Fe is slightly larger than the 
one of O, because our adopted prescription for the SNIa 
rate produces a smaller O /Fe ratio in the inner disk than 
in the outer one. The gradient of secondary N is always 
steeper than the one of O, but since we do not include N 
production from intermediate mass stars or WR stars in 
this calculation, this result serves merely for illustration 
purposes. 

4. Comparison to observations 

In this section we compare our results to the observed 
abundance profiles of various elements across the galactic 
disk. Tracers of abundance profiles include emission-line 
objects (HII regions and planetary nebulae), as well as 
stars and stellar associations (B-stars and open clusters, 
respectively). In most cases, those tracers are relatively 
young objects (HII regions, B-stars and planetary nebu- 
lae of type I), younger than ~ 1 Gyr; they provide then 
information about the current abundance profile of the 
Milky Way disk. In other cases (planetary nebulae of type 
II or III, open clusters), the objects involved are several 
Gyr old and provide information about the past status of 
the disk (albeit with much larger uncertainties than in the 
former case). 

In Sect. 4.1 we present briefly the available observa- 
tional data from various sources. In Sect. 4.2 we compare 
the data for young objects to our results at T = 13.5 Gyr; 
we show, in particular, how the results for C (and to a 
much lesser extent, N and O) may be affected by the M92 
metallicity-dependent yields of massive, mass losing stars. 
In Sect. 4.3 we discuss the profiles of the corresponding 
abundance ratios of elements to Oxygen. In Sect. 4.4 we 
discuss the history of the abundance profiles, comparing 
our results to observations of old objects (of rather uncer- 
tain ages). 

4-1. Observational data 

One of the first comprehensive optical surveys of HII re- 
gion abundances was performed by Shaver et al. (1983). 
They have found strong abundance gradients for N/H, 
O/H and Ar/H in the range of galactocentric distances 
Rg = 5—12 kpc. Those results were later confirmed by 
the surveys of Fich & Silkey(1991) and Vilchez & Este- 
ban (1996), extending the data up to distances of Rg ~ 

17 kpc. Those authors suggested that the N/H and O/H 
gradients show a tendency for flatenning between 11 and 

18 kpc. Taking advantage of the small extinction in the in- 
frared, Simpson et al. (1995) and Afflerbach et al. (1997) 
studied objects towards the center of the Galaxy. Simp- 



son et al. (1995) observed 12 HII regions in Rg = 0—10 
kpc and found a somewhat better fit with a step function 
than with a smoothly decreasing one. On the contrary, 
after adding five outer Galaxy HII regions (Rg = 13 — 
17 kpc), Rudolph et al. (1997) concluded that the single 
slope profiles for N/H and S/H are more likely than step 
functions. The existence of a gradient in the oxygen abun- 
dance profile obtained from HII regions is not in doubt at 
present, but its magnitude has been recently challenged by 
the work of Deharveng et al. (2000), who report a value 
about half as large as previous measurements. 

Observations of B-type stars in young clusters and as- 
sociations in the late 90ies confirmed the existence of an 
abundance gradient. Smartt & Rollcston (1997) found an 
oxygen abundance gradient similar to the one obtained 
in HII regions (except for the work of Deharveng et al. 
2000), in the galactocentric distance range Rg = 6 — 18 
kpc. This was confirmed by Gummersbach et al. (1998) 
who detected significant abundance gradients for C, O, 
Mg, Al and Si in galactocentric distances Rg = 5 — 14 
kpc. 

Several works have been devoted to the study of the 
abundance patterns of planetary nebulae (PN, see e.g. Ma- 
ciel and Quiroza 1999 and references therein). According 
to a classification scheme originally suggested in Peimbert 
(1978) and revised in Pasquali and Perinotto (1993), plan- 
etary nebulae are subdivided into types I, II and III. PNI 
have He/H > 0.125 and log(N/0) > -0.3, arc located in 
the galactic thin disk and their progenitors arc probably 
associated with stars in the mass range 2.5 — 8 M Q . PNII 
belong also to the thin disk, are not particularly enriched 
in nitrogen or He and are thought to evolve from stars 
in the 1.2 — 2.5 M Q range. Finally, PNIII are associated 
with the galactic thick disk and thought to originate from 
low mass progenitors, in the 1 — 1.2 M Q range. The esti- 
mated progenitor masses imply ages of <1 Gyr for PNI, 
1 - 8 Gyr for PNII and >8 Gyr for PNIII, respectively. 
These differences should allow, in principle, to use PN as 
tracers of the chemical evolution of the galactic disk. How- 
ever, the mass and age differences between the various PN 
types are not quite well defined, as we shall sec in Sect. 
4.3. Moreover, all PN are expected to be auto-enriched in 
products of the CN cycle (i.e. N-rich and C-poor), while 
PNI (especially those originating from relatively massive 
progenitors, in the 5 — 8 M Q range) are probably auto- 
enriched in products of the NO cycle (N-rich and O-poor) . 
For those reasons, we shall not consider at all in the fol- 
lowing the abundances of C and N from PN of all types 
and the abundances of O from PNI. For the purpose of 
this work, we shall use PNI (along with HII regions and 
B-stars) as tracers of the young disk in Sects. 4.2 and 4.3. 
PNII and PNIII will be discussed in Sect. 4.4. 

Table 1 summarizes the currently available observa- 
tional data on the abundance profiles across the Milky 
Way disk. Two points should be noticed: 
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1) Contrary to other authors, we did not consider data 
on Fe from open clusters, since the situation is rather 
uncertain at present: for instance, in the solar neighbor- 
hood, open cluster data show no evolution of the Fe abun- 
dance with age (e.g. Friel 1999), contrary to the familiar 
age-metallicity relation suggested by data for F-stars (Ed- 
vardsson et al. 1993). In view of that dicrepancy, we chose 
not to consider Fe in this work (although we do calculate 
its abundance profile, as shown in Fig. 4). 

2) Data in Table 1 are adopted directly from the ref- 
erences listed and no attempt has been made to homoge- 
nize them by recalculating the abundances in a consistent 
way. In general, differences in techniques and atomic data 
employed produce an abundance scatter which is smaller 
than observational uncertainties in the line strengths. 
Therefore, we believe that direct comparison between our 
models and those inhomogeneous data can still provide 
statistically meaningful results. 

4-2. Current disk abundance profiles 

The observational data of Table 1 concerning "young" ob- 
jects (HII regions, B-stars, PNI) are plotted in Fig. 6, 
along with our model results, obtained with the WW95 
yields, at T = 13.5 Gyr. In Fig. 7 we also plot the same 
data for He, C, N and O and we show the correspond- 
ing results obtained with both WW95 and M92 yields; 
as discussed in Sect. 2.2, the metallicity dependenence of 
massive star winds has an effect on the yields of those 
elements. 

Since the prescriptions for the radial dependence for 
infall and SFR of the BP99 model adopted here were such 
as to reproduce the observed oxygen abundance profile, 
we start the description of our results with this element. 

Oxygen: As can be seen from Table 1, data from both 
HII regions and B-stars suggest an abundance gradient 
dlog(0/H)/dR ~ —0.07 dex kpc -1 (see the recent review 
by Smarrt 2000 on B-stars). Maciel and Quiroza (1999) 
include PNII (despite the fact that these objects are, in 
principle, older than 1 Gyr) and find an average gradient 

of dlog(0/H)/dR 0.065 dex kpc^for all tracers of the 

"young" population. However, Deharvenget al. (2000), af- 
ter a consistent analysis of their own data on HII regions, 
as well as of those of previous works, conclude that the 
gradient should be ~ 40% smaller than generally thought, 
i.e. -0.04 dex kpc" 1 . The data of Deharveng et al. (2000) 
concern the galactocentric distance Rq = 5—15 kpc, i.e. 
the same as the one of the work of Shaver et al. (1983) 
on HII regions or studies on B-stars; the difference in the 
resulting abundance gradient cannot then be attributed in 
the studied galactocentric distance range. As stressed by 
Deharveng et al. (2000), their derived O/H abundances de- 
pend strongly on their two-temperature HII region model 
(one temperature for the high excitation ++ zone and 
another temperature for the low excitation + zone); an 



alternative HII region model (with temperature fluctua- 
tions) would lead to larger O/H abundances. 

As explained in Sect. 2, our model is based on the ra- 
dial variation of the SFR (Eq. 3) and infall rate (Eq. 1). 
A different prescription would lead to different results for 
the abundance gradients (as in e.g. Prantzos and Aubert 
1995, who used a radially independent infall time-scale 
and found a smaller gradient of dlog(0/H)/dR ~ —0.03 
dex kpc" 1 ). Also, by adopting the same prescriptions, but 
assuming radial inflows, one would obtain smaller values 
for the abundance gradient (e.g. Tsujimoto et al. 1995, 
Portinari and Chiosi 2000). Obviously, the magnitude of 
the abundance gradient is crucial in fixing the parame- 
ters entering the current phcnomcnological chemical evo- 
lution models. In particular, a strong gradient does not 
support the idea of important radial inflows induced by a 
galactic bar (such as those found e.g. in the calculations 
of Fricdli and Benz 1995); the Milky Way bar must then 
have formed relatively recently and/or played a negligible 
role in driving gaseous flows in the disk. 

The only published successful chemodynamical model 
for the Milky Way (Samland et al. 1997) predicts an 
oxygen abundance gradient that flattens considerably in 
the outer disk. Despite some claims for such a flatenning 
(Vilchez and Esteban 1996), no convincing observational 
evidence exists up to now (see the discussion in Deharveng 
et al. 2000). Our model predicts no flatenning of the gra- 
dient in the outer disk. It does predict a small flatenning 
in the inner disk. As explained in Prantzos and Aubert 
(1995), this is due to the fact that all metal abundances 
are diluted in the inner disk by the late ejection of metal- 
free material by the numerous low-mass, long-lived stars 
of the first stellar generations; the ejection rate of that 
material at late times is larger than the metal production 
rate, since most of the gas in the inner disk has been con- 
sumed at that time. This is not the case in the outer disk, 
which is formed late. 

As can be seen in Fig. 7, the M92 yields lead to results 
for oxygen that are almost indistinguishable from those of 
WW95. Only in the inner (more metal-rich) zones a slight 
difference is obtained in the oxygen abundance profile. As 
already discussed in Prantzos et al. (1994) the metallicity 
dependent stellar winds have a negligible impact on the 
oxygen yields. 

Finally, another aspect of the observed oxygen abun- 
dances (Fig. 6) has been discussed in several places: all 
"young" objects in the local disk have lower oxygen abun- 
dances (in the range O/H ~ 3.2—5 10~ 4 ) than the Sun 
((O/H) = 6.8to.g 10~ 4 , Grevesse & Sauval 1998). For 
instance, Gies and Lambert (1992) derive an oxygen abun- 
dace 0/H~ 4.8 10~ 4 for local B stars, while Esteban et 
al.(1998) give O/H = 5.2^ 10~ 4 for Orion nebula. This 
is difficult to understand in the framework of conventional 
models of chemical evolution, where metallicity increases 
monotonically with time. The idea that the Sun was born 
in the metal-rich inner disk and subsequently migrated 
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Fig. 6. Present day abundance profiles in the Milky Way disk (obtained in the framework of the model presented in 
Sect. 2 and with the WW95 metallicity dependent yields) and comparison with observations. Data sources are given in 
Table 1. The HII region data of Deharveng et al. (2000) suggest an oxygen abundance gradient by 40% less steep than 
the commonly accepted one of dlog(0/H)/dR ~ —0.065 dex/kpc. The typical observational uncertainties for different 
elements and data sets are shown with error bars in the bottom left corner of each panel. 



outwards (Wielen et al. 1996), has been recently rejected 
by Binney and Sellwood (2000), on the grounds of dynam- 
ical arguments. Thus, at present, there is no satisfactory 
explanation for the "super-metallicity" of the Sun relative 
to young objects in the solar neighborhood. 

Carbon: There are very few studies of carbon abun- 
dances, either in HII regions or B-stars. Based on only two 
nearby HII regions (Orion and M17), Peimbert et al.(1992) 
first derived an abundance gradient of dlog(C/H)/dR ~ 
— 0.08±0.02 dex kpc~Mn the solar vicinity. Adding an- 
other HII region (M8) to that sample, Esteban et al. 
(1999) obtain -0.133±0.022 dex kpc _1 . However, both 
those studies concern a very limited range of galactocen- 
tric distances (6 — 9 kpc), and the derived abundance 
gradients cannot be considered as representative of the 
disk as a whole. On the other hand, the study of Gum- 
mersbach et al. (1998) suggests a rather small C gradient 
for B-stars in the range Rg = 6 — 12 kpc: dlog(C/H)/dR 



~ —0.045 dex kpc -1 ; this is half as steep as the oxygen 
abundance gradient derived for those same B-stars in that 
study. If the Gummersbach et al. (1998) value is correct, it 
should then be difficult to understand the corresponding 
disk profile of the C/O ratio, which should decrease in the 
inner disk. Indeed, observations of low mass stars in the 
solar neighborhood show that C/O increases with metal- 
licity (e.g. Gustafsson et al. 1999 and references therein), 
and we would expect to see the same effect in the inner, 
metal-rich disk. Hibbins et al. (1998) performed a differ- 
ential analysis of C and N abundances of B-stars in the 
outer disk (10 — 17 kpc). They found that C correlates 
with N, but not with O. However, their absolute values of 
C/H for nearby stars are systematically lower by a factor 
2^3 than those of other studies; they suggest then that 
their data (and data coming from different samples, in 
general) should be used with caution in studies of galactic 
chemical evolution. 
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Fig. 7. Abundance profiles of He, C, N and O, at T = 
9 Gyr (dashed curves) and T = 13.5 Gyr (solid curves), 
obtained with the WW95 yields (thin curves, same as e.g. 
in Fig. 3), and the M92 yields (thick curves), in the frame- 
work of the model presented in Sect. 2. It can be seen that, 
compared to the WW95 yields, the M92 yields lead to con- 
siderably different results for C, to slightly different results 
for N, and to quasi- identical results for He and O. In par- 
ticular, the M92 yields may explain completely the carbon 
abundance profile, with no need for any contribution from 
intermediate mass stars; on the contrary, nitrogen obser- 
vations suggest that another source is mandatory. 



In summary, the carbon abundance profile in the Milky 
Way disk is poorly determined. The few available data 
sets cover limited galactocentric distance ranges and are 
not derived in a consistent way. If all the data are plotted, 
as in Fig. 6 and 7, a carbon abundance gradient larger 
than the one of oxygen appears; but this may well be an 
artifact. We note that in his recent review, Smartt (2000) 
suggests a gradient of dlog(C/H)/dR ~ —0.07 dex kpc -1 , 



based on unpublished work of Rolleston ct al. (2000) on 
B-stars. 

Our model results (Fig. 6) show a carbon abundance 
gradient of dlog(C/H)/dR ~ —0.06 dex kpc -1 , i.e. com- 
parable to the one of oxygen, when the WW95 yields 
are used. However, the absolute value of the local carbon 
abundance at the Sun's birth is underproduced in that 
case (by a factor of ~2); another C source is then required. 
Intermediate mass stars (IMS) are well known net produc- 
ers of carbon, but in recent years, several works (Prantzos 
et al. 1994, Carigi 1994, Gustaffson et al. 1999, Henry et 
al. 2000) suggested that the M92 metallicity dependent 
yields of massive stars fit better the observed C/O abun- 
dance patterns in extragalactic HII regions and low-mass 
stars in the solar neighborhood. This is confirmed in Fig. 7, 
where it can be seen that the use of the M92 yields allows 
indeed to reproduce the absolute abundance of carbon in 
the solar neighborhood, with no need for a contribution by 
IMS. In that case, the resulting C abundance gradient is 
much steeper: —0.086 dex kpc -1 . Clearly, the abundance 
profiles of carbon vs. oxygen in the Milky Way disk are 
crucial in any attempt to evaluate the role of massive star 
winds in the production of carbon. 

Nitrogen: There is a large number of works concern- 
ing N abundances in HII regions; the N abundance profile 
seems to be steeper, in general, than the corresponding O 
profile (see Table 1). Observations of B-stars (Gummers- 
bach ct al. 1998) support this conclusion, but the error 
bars are much larger in that case. Smarrt (2000) suggests 
a gradient of dlog(N/H)/dR 0.08±0.002, i.e. compat- 
ible with the one of oxygen. Our models also produce a 
steeper gradient of N (with respect to that of O), with the 
yields of both WW95 (Fig. 6) and M92 (Fig. 7). As can 
be seen in Fig. 7, the differences resulting from the use 
of those two sets of yields are rather small (slightly larger 
than in the case of oxygen, but considerably smaller than 
in the case of carbon). Two important points should be 
noticed: 

i) The absolute value of the N abundance obtained 
with both sets of yields is too low compared with observa- 
tions in the solar neighborhood and in the disk. Contrary 
to the case of carbon, the M92 yields are not sufficient to 
account for the current abundance of N in the Milky Way. 
Unless the N yields are seriously underestimated in both 
WW95 and M92 (which is improbable, since they con- 
sist of the sum of the initial C+N+O), another N source 
is required. Intermediate mass stars are the obvious can- 
didate, but their yields are notoriously difficult to esti- 
mate (in view of the many uncertainties concerning mass 
loss rates, "hot-bottom burning", etc.). Our calculations 
clearly show what is the magnitude of the expected con- 
tribution of IMS, in order to complement the (presum- 
ably better understood) yields of massive stars: on aver- 
age, IMS have to produce 3-4 times more N than massive 
stars. 
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ii) We stress again (see also Sect. 3) that the flatenning 
of the N abundance profile obtained in our calculations 
for the inner disk is due to the fact that, since yields of 
WW95 and M92 are available only up to Z = Z©, for 
higher metallicities we use the yields at Z = Z©; thus, in 
our models N is produced in the inner disk (where high 
metallicities are developed) with the same yield always, 
i.e. as a primary element and its abundance follows the 
one of oxygen. Clearly, this is an artifact of the calculation, 
due to the lack of appropriate input data; if N were treated 
correctly, as a secondary, its abundance profile would not 
flatten in the inner disk [Notice: we could simply scale the 
N yields of WW95 with metallicity, but it would be more 
difficult to do the same thing with the M92 yields, since 
the effect of the stellar wind on the yield is not simply 
proportional to metallicity] . 

Magnesium, Aluminium, Silicon: The abundance 
profiles of those elements have been studied by Gummesr- 
bach ct al. (1998), who observed B-stars in galactocentric 
distances 5 — 14 kpc. In view of the large uncertainties 
(see Table 1), all those gradients can be considered as 
compatible with the oxygen gradient of —0.07 dex kpc -1 , 
at the 1 a level; this is particularly true for Mg. On the 
other hand, taken at face value, the Al profile seems to be 
considerably steeper than the one of Si (—0.045 dex kpc -1 
vs. —0.107 dex kpc" 1 , respectively). This is rather surpris- 
ing, taking into account that Al is an odd-Z element and 
its yield depends slightly on metallicity (the "odd-even" 
effect, see Fig. 1). Because of this effect, our models pro- 
duce a steeper gradient for Al than for Si. The latter is 
comparable to the one of oxygen, as expected. 

We note that in his recent review Smarrt (2000) 
suggests gradients of similar magnitude for Al and Si 
(-0.05±0.02 dex kpc" 1 and -0.06±0.01 dex kpc" 1 , re- 
spectively), based on B-star data. In that case the prob- 
lem is alleviated, although not completely solved. We also 
note that low-mass stars in the solar neighborhood do not 
exhibit the theoretically expected behaviour of the Al/O 
abundance pattern (e.g. Goswami and Prantzos 2000 and 
references therein) . It may well be then that the odd-even 
effect has been overestimated in the WW95 yields of Al. 

In summary, the observed abundance profiles of Al and 
Si run opposite to theoretical expectations, but in view of 
their large uncertainties, it is difficult to draw firm con- 
clusions on stellar nucleosynthesis. 

Neon, Sulphur, Argon: These elements have been 
observed through their emission lines in both HII regions 
and PN of all types. In general, HII region abundances 
suggest a radial profile similar to the one of oxygen for 
all three elements, i.e. a gradient of —0.07 dex kpc -1 (the 
flat S profile of Shaver et al. 1983 is an exception to this 
general agreement). Data from PNI and PNII support gra- 
dients of this magnitude also for S and Ar, but suggest a 
smaller gradient for Ne (—0.036 dex kpc -1 , according to 
Maciel and Quiroza 1999). 
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Fig. 8. Present day abundance gradients in the Milky 
Way disk: summary of models vs. observations. Model re- 
sults are shown by open symbols (asterisks, for the WW95 
yields, connected by a solid curve, and pentagons for the 
M92 yields, only for C, N and O). Observations (Table 1 
for references) are shown by filled symbols. Notice the large 
discrepancy between the carbon gradients observed in B- 
stars and HII regions, respectively. An average gradient 
of — 0.07±0.02 dex kpc _1 is compatible with most of the 
data and the model results. Most of the observed abun- 
dance gradients are rather well reproduced by the model, 
but the case of Al is problematic (see text). 



Our model leads to S and Ar abundance gradients sim- 
ilar to the one of O (~ —0.05 dex kpc -1 in both cases). 
In the case of Ne, the (unexepected) small metallicity de- 
pendence of the WW95 yields (Fig. 1) leads to a steeper 
profile: dlog(Ne/H)/dR ~ —0.08 dex kpc -1 , exactly as in 
the case of Mg. 

In summary, the observed S and Ar abundance profiles 
are in agreement with theoretical expectations, while Ne 
remains problematic, both theoretically and observation- 
ally. 

Helium: For the sake of completeness, we present here 
the He profiles of our models. As can be seen in Fig. 7, 
there is no difference between the profiles obtained with 
the WW95 and M92 yields. In both cases a flat He/H 
profile is obtained. However, the lack of IMS yields in our 
model does not allow to draw conclusions about He. 

The observational situation concerning the He abun- 
dance profile is not clear. The flat profile of the He + /H + 
ratio found by Shaver et al. (1983) was confirmed by the 
recent work of Deharveng et al. (2000). However, the true 
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Hc/H gradient depends also on the unknown amount of 
the corresponding neutral species. The data on Fig. 6 
and 7 are from the work of Gummersbach et al. (1998) 
on B-stars. The error bars are quite large (±0.3 dex for 
log(Hc/H)) and do not allow to draw any conclusion on 
the existence of an abundance gradient for He. 

Fig. 8 summarizes the discussion of this section, con- 
cerning observed and calculated abundance profiles in the 
Milky Way disk. If the values of Smartt (2000) are adopted 
for the C and Si gradients from B-stars, then there is sat- 
isfactory agreement between the model and observations 
for all elements but Al (taking into account error bars); 
He and N require another source to account for their ab- 
solute abundances, even if the gradient has the correct 
value. But the main "message" of Fig. 8 is that homoge- 
neous data sets are required for a meaningful comparison 
between theory and observations. 

4-3. Abundance ratios along the Galactic disk 

Abundance ratios between metals are, in principle, more 
reliable tracers of the chemical evolution than absolute 
abundances, since they do not depend on the star forma- 
tion efficiency and they allow to identify if e.g. an element 
has a secondary origin, whether it is produced in long-lived 
sources (low-mass stars or SNIa) or whether it is affected 
by the "odd-even" effect, etc. These properties have been 
widely used as diagnostic tools of the history of the halo 
+ solar neighborhood system (se e.g. Pagel 1997). In prac- 
tice, the situation is complicated due to various theoretical 
and observational uncertainties (e.g. Goswami and Prant- 
zos 2000 and references therein), except for a few trivial 
cases. 

This is also the case for the Milky Way disk, as can 
be seen in Fig. 9: of the nine element abundance ratios 
to oxygen displayed in the figure, not a single one shows 
any clear trend with galactocentric distance. Although in 
most cases this is expected on theoretical grounds (e.g. for 
Ne, Mg, Si, S, Ar, which are primary elements), for oth- 
ers (N and, to a smaller extent, He and Al) this is rather 
surprising. Even worse is the large scatter obtained for all 
elements and at all galactocentric distances. There is a 
striking difference with F stars in the solar neighborhood, 
which display very small dispersion in their abundance ra- 
tios (Edvardsson et al. 1993). The inhomogeneity of the 
data sets displayed in Fig. 9 contributes certainly to that 
scatter. However, even in cases where only one tracer is in- 
volved (B-stars for Mg, Al and Si), the dispersion is rather 
large: it is about twice as large as the typical uncertainty 
of individual abundance ratios, estimated here to be ±0.2 
dex (see also Henry and Worthey 1999). 

A meaningful comparison between theory and observa- 
tions is difficult in such conditions. However, Fig. 9 allows 
to draw some conclusions: 

i) The theoretical abundance ratios of Si/O, S/O and 
Ar/O are ~ constant across the disk, as expected for pri- 



maries; their absolute value is always ~ solar and compat- 
ible with available observations. The WW95 yields repro- 
duce well the current average abundance ratios of those 
elements. 

ii) The theoretically expected "odd-even" effect for 
Al/O is not manifested in the available observational data, 
i.e. no trend of Al/O with galactocentric distance (or 
metallicity) is observed. The WW95 yields reproduce cor- 
rectly the current local Al/O ratio (~ solar), and lead to 
a small dependence on metallicity (or radius), because of 
the "odd-even" effect. 

iii) The WW95 yields underproduce the observed 
Ne/O and Mg/O ratios (as anticipated from Fig. 1); this 
is also the case for the solar neighborhood, as discussed 
in Goswami and Prantzos (2000). The discrepancy is not 
large, taken into account the various uncertainties of stel- 
lar nucleosynthesis calculations (see Prantzos 2000 for a 
review) as well as in observational data. Our Ne/O data 
come mainly from PN I, where this ratio is usually deter- 
mined under the assumption that Ne/O = Ne ++ /0 ++ . 
The Ne ++ /0 ++ values probably are upper limits to the 
Ne/O value for central stars hotter than 50,000 degrees, 
due to the presence of charge exchange reaction ++ + H° 
=>■ + + H+ that allows some + to coexist with Ne ++ ; 
this is particularly the case for PNs of type I with low 
density (Pcimbert et al. 1995). Therefore the best values 
from type I PNs are given by the lower envelope of the 
data presented in Fig. 9 (M. Pcimbert, private communi- 
cation). On the other hand, the WW95 yields of Ne and 
Mg manifest an unexplained dependence on metallicity. 
Put together, these discrepancies point to some problems 
in the WW95 yields of Mg and Al. One possibility is that 
the extent of the C-shell (where both elements are mainly 
synthesized) is underestimated by the Ledoux criterion for 
convection employed in the WW95 calculations. 

iv) The WW95 yields underproduce the local C/O ra- 
tio, while the M92 yields of massive stars reproduce it 
fairly well (with no need for a C contribution by IMS). 
However, if WR stars are the main producers of C, one 
expects this ratio to decrease with galactocentric distance 
(as indicated by the thick curve for C/O in Fig. 9), and 
this trend is not seen in currently available data. On the 
other hand, if IMS stars produce the bulk of C as primary, 
no variation of the C/O ratio with galactocentric distance 
is expected. It should be noted that gradients of C/O (in 
fact, variations of C/O with O/H) have been observed in 
extragalactic HII regions (e.g. Garnett et al. 1999) and are 
also expected in the case of the Milky Way disk. Thus, ac- 
curate determination of the C/O ratio across the Milky 
Way disk is crucial in determining the roles of IMS and 
Wolf-Rayet stars in carbon production. [Note: the flatten- 
ing of the theoretical C/O profile in the inner disk is due to 
our use of M92 yields for Z = Z Q even when the metallic- 
ity is larger than solar; in principle, the C/O ratio should 
continue increasing in the inner disk. C and O yields for 



J.L. Hou, N. Prantzos and S. Boissier: Abundance gradients and their evolution in the Milky Way disk 13 




R G (kpc) 

Fig. 9. Profiles of abundance ratios versus oxygen (normalised to the corresponding solar ratio: 
[X/0]=log((X/0)/(X/0)©)) in the Milky Way disk. Observations are shown and compared to model results of this 
work. In all panels, the thin curves indicate results obtained with the WW95 yields, while the thick curves on the left 
side panels indicate results obtained with the M92 yields (leading to larger He/O, C/O and N/O abundance ratios). 
[Note: The flattening of theoretical profiles of abundance ratios in the inner disk is due to our use of yields for Z=Z Q 
even at metallicities larger than solar, because yields for such metallicities are not available in WW95 and M92]. Data 
(from references given in Table 1) concern HII regions (triangles), PNI (crosses) and B-stars (pentagons). The typical 
observational uncertainties for different element ratios are shown with error bars in the upper left corner of each panel. 
Notice that for Ne/O, the best values from type I PNs are given by the lower envelope of the data (see discussion in 
Sect. 4.3). 



metallicities higher than solar are required for a proper 
calculation] . 

v) Both WW95 and M92 yields underproduce the N/O 
ratio in the disk and lead to a N/O profile declining with 
radius. IMS stars are expected to be the main N producers 
in the disk. The dispersion of presently available data on 
N/O in the disk does not allow to conclude whether they 
produce N as primary (i.e. through hot-bottom burning), 
or as a secondary. As with the case of C, observations of 
N/O in extragalactic HII regions show that N behaves as 
secondary at high metallicities (Henry et al. 2000). This 
does not seem to be the case in the inner regions of the 
Milky Way disk, at least with currently available data. 



vi) The observed He/O ratio in B-stars is higher than 
the corresponding solar value at all galactocentric dis- 
tances. This is obviously related to the fact that the O 
abundances of those young objects are (somewhat surpris- 
ingly) systematically lower than solar. Since no satisfac- 
tory explanation for that exists up to now (at least in the 
framework of conventional chemical evolution models), we 
do not expect our results to match the average observed 
value of Hc/O. On the other hand, the B-star data show 
a flat He/O profile (as a result of the similar He/H and 
O/H gradients shown on Fig. 6 and 7). In our calculations, 
the amount of He produced by massive stars alone (either 
with the WW95 or the M92 yields) is negligible relative 
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Fig. 10. Abundance profiles in the Milky Way disk at various epochs according to our model and comparison to 
objects of (hopefully) corresponding ages. From left to right, profiles are shown for young (column A), intermediate 
(column B) and old age (column C) objects. In Case A belong B-stars, HII regions and planetary nebulae of type I; 
in Case B belong planetary nebulae of type II; and in Case C belong planetary nebulae of type III (see Sect. 4.4 for a 
discussion on ages). The corresponding model results are shown at time T = 12.5—13.5 Gyr or age = 0—1 Gyr (shaded 
area in A), T = 5.5—12 Gyr or age = 1.5—8 Gyr (shaded area in B) and T = 0.5—5.5 Gyr or age = 8—13 Gyr (shaded 
area in C). References for data are given in Table 1. 



to the primordial one of the infalling gas. Thus, the final 
He/H profile is flat (Fig. 6 and 7) and the final He/O pro- 
file (Fig. 9) is decreasing in the inner disk; obviously, the 
contribution of intermediate mass stars is mandatory to 
account for He observations. 

4-4- Evolution of abundance gradients 

As already mentioned in Sect. 4.1, abundances in plan- 
etary nebulae of various types (I, II and III) allow, in 
principle, to follow the history of the abundance profile 
of the Milky Way disk. However, the progenitor masses 
and lifetimes of PNs are not well known. For instance, 
Allen et al. (1998) assume that both PNII and PNIII have 
progenitors of similar mass range (M<2 M©) and classify 
them in terms of their kinematical properties: PNII have 



peculiar velocities vp < 60 km/s and PNIII have vp > 
60 km/s. To account for these properties as well as for 
the relative frequency of the various PN types, Allen et 
al. (1998) propose a "dynamical" model for the Milky 
Way disk, simulating orbital diffusion. They find good 
agreement with observed gradients of PN (classified ac- 
cording to that scheme), provided that not all IMS go 
through the PN stage. In that same paper, they also ex- 
plore the "conventional" scheme, classifying PN in terms 
of progenitor mass: 1.3<M/M©<8.4 and r < 3 Gyr for 
PNI, 0.9<M/M©<1.3 and 3< r < 9 Gyr for PNII and 
0.8<M/M©<0.9 and t > 9 Gyr for PNIII. They find that 
their model fails to reproduce observations of the various 
PN types in that case. 

In our model we have no information on the kinemat- 
ical properties of the gas or the stars. We adopt then the 
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classification scheme of Pasquali and Perinoto (1993) for 
PN, already presented in Sect. 4.1. We are aware that this 
scheme is not necessarily the most appropriate one (uncer- 
tainties in ages may be larger than stated and kinematics 
may play a non negligible role), but we adopt it for the 
sake of comparison to our model results. This comparison 
is presented in Fig. 10. It can be seen that: 

i) The model describes relatively well the O, Ne, S 
and Ar abundance gradients of young objects (left hand 
panels, column A), as already discussed in the previous 
sections. However, the observed scatter at all galactocen- 
tric distances is much larger than the range of values given 
by the model for the last 1 Gyr (the shaded region in Fig. 
10). It should be noted that the observed scatter (<~ 0.6 
dex) is larger than the corresponding one for young stars 
in the solar neighborhood (e.g. Garnett and Kobulnicky 
2000). Probably, systematic effects and uncertainties in 
distance estimates contribute largely to the observed scat- 
ter of "young" tracers in Fig. 10. 

ii) Model predictions are compatible with observations 
for O and Ar in PNII, as can be seen in the middle panels 
of Fig. 10 (column B). In fact, there is very little difference 
between the data for PNII and those for PNI. Predictions 
for S are slightly above the data points, while predictions 
for Ne are clearly below the corresponding data. The rea- 
son for the latter discrepancy is obviously the unexplained 
underproduction of Ne in the low mctallicity yields of 
WW95, already discussed in Sects. 2.2 and 4.3. In the 
other 3 cases (O, S and Ar) there is fairly good agreement 
between theory and observations concerning all proper- 
ties of the observed abundance profiles (absolute values, 
gradient, scatter). 

iii) Observations show that, on average, abundances 
in PNIII are lower and present a larger dispersion than in 
PNII. Both features are relatively well reproduced by our 
models, as can be seen on the right hand panels in Fig. 
10 (column C), with the exception of Ne (for the reasons 
already mentioned in the previous paragraph). 

Our models, as some other models of this kind (Molla 
et al. 1997, Portinari and Chiosi 1999) suggest that abun- 
dance gradients are steeper for older objects. Other works 
(e.g. Chiappini et al. 1999) reach an opposite conclusion, 
although they are based on the same assumption, namely 
an "inside-out" formation of the disk; clearly, the adopted 
time-scales for star formation and infall play an important 
role in this discrepancy. Although our model results are 
compatible with available data on PN, the large scatter 
in those data does not allow to conclude on the temporal 
variation of the gradients (see Maciel and Quireza 1999). 
However, our model makes a new, and perhaps testable, 
prediction: if the observed abundance scatter at a given 
galactocentric distance is to be attributed, at least par- 
tially, to intrinsic age differences between the tracers (e.g. 
PN) then we expect that scatter to be smaller in the inner 
disk than in the outer one (as indicated by the extent of 
the shaded areas in Fig. 10). 



5. Summary 

Based on the successful chemical evolution model for the 
Milky Way disk developped by BP99, we calculated the 
corresponding abundance profile of elements up to the Fe- 
pcak. For that purpose we used the metallicity dependent 
yields of WW95 for massive stars evolving at constant 
mass, along with those of Iwamoto et al. (1999) for SNIa. 
We also explored the yields of M92 for massive mass losing 
stars, concerning the elements He, C, N and O. We delib- 
erately ignored yields from intermediate mass stars (IMS), 
our purpose being to check to what extent such stars are 
indeed required to explain observed abundance patterns. 
However, we did take into account ejecta from IMS, as- 
suming that their net yield is zero for all elements; as ex- 
plained in Goswami and Prantzos (2000), this procedure 
is crucial to ensure that the absolute metal abundances, 
expressed as X/H, are correctly evaluated. 

We compiled a large sample of observational data, con- 
cerning abundance profiles of He, C, N, O, Ne, Mg, Al, 
Si, S and Ar in the Milky Way disk (Table 1). Most of 
the tracers are young objects (B-stars, HII regions, PNI), 
while PNII and PNIII abundances trace earlier stages of 
the Galaxy evolution (but with considerable uncertainties 
in the corresponding ages). Contrary to other authors, we 
did not consider data on Fe from open clusters, since the 
situation is rather uncertain at present (see Sect. 4.1). We 
simply note that our prescription for the rate of SNIa (ma- 
jor Fe producers), leads to a Fe abundance profile slightly 
steeper than the one of O in the disk (Fig. 4 and 5). 

The main results of the model may be summarized as 
follows: 

(i) We obtain abundance gradients for all elements 
from He to Zn. For primary elements, like O, the the- 
oretical value of the abundance gradient at the present 

epoch (T = 13.5 Gyr) is: dlog(X/H)/dR 0.06 dex 

kpc -1 . This value is compatible with most observational 
data concerning young objects in the disk (see Fig. 8 and 
Table 1). We note, however, the "puzzling" conclusion of 
Dcharveng et al. (2000), pointing to an oxygen gradient 
smaller by about 40% than the commonly accepted value; 
if their conclusion is confirmed, our model parameters (the 
ratio of star formation to infall timescale as a function of 
galactocentric radius) should be appropriately modified. 
For secondary or "odd-Z" elements (like N and Al, re- 
spectively) we obtain slightly larger gradients; this is also 
the case for C when the M92 yields arc used. Such values 
are marginally compatible with available observations. We 
find that the M92 yields can account for the totality of 
carbon production across the disk with no need for a con- 
tribution by IMS, in agreement with other recent studies 
(Prantzos et al. 1994, Carigi 1994, Gustafsson et al. 1999, 
Henry et al. 2000a, 2000b). On the contrary, the majority 
of N is produced by another source, most probably IMS. 

(ii) Current observations show no trend of the abun- 
dance ratio X/O with galactocentric distance (or metal- 
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licity), for any element X. This is unexpected in the case 
of N and C, which do show such a trend in extragalac- 
tic HII regions where their abundance ratio to oxygen in- 
creases with metallicity (e.g. Henry and Worthey 1999). 
Our model shows a flat profile of S/O, Si/O and Ar/O (in 
agreement with observations) and a slowly decreasing ra- 
tio of Ne/O, Mg/O and Al/O with galactocentric distance 
(due to a "problematic" or overestimated metallicity de- 
pendence of the WW95 yields for those elements). They 
also show that C/O should decrease with galactocentric 
distance in the case of M92 yields (but this should not be 
the case if C is mainly produced as a primary by IMS). 
Clearly, precise observations of the C/O ratio across the 
disk are required in order to decide about its main produc- 
tion site. Similar conclusions hold for N: the M92 yields 
cannot account for its abundance profile, but is not clear 
whether primary or secondary production in IMS is the 
dominant mechanism. Observations of extragalactic HII 
regions suggest that the latter mechanism dominates in 
high metallicity regions (Henry et al. 2000b); if this is 
true, then an important N/O gradient should be observed 
in the Milky Way disk. 

(iii) The evolution of abundance gradients provides a 
strong constraint for Galactic chemical evolution models. 
Despite a considerable amount of observational and the- 
oretical work (e.g. Henry and Worthey 1999; Tosi 2000 
and references therein) , the question has not yet been set- 
tled. Observations of planetary nebulae (PN) of different 
ages could, in principle, help in that respect (Pasquali and 
Perinoto 1993, Maciel and Quireza 1999), but the ages 
and distances attributed to PN are not sufficiently well 
known at present. Our model predicts a steady flattening 
of the gradients with time, due to the adopted "inside- 
out" formation scheme for the disk. Such an evolution is 
also found in other works (e.g. Molla et al. 1997, Portinari 
and Chiosi 1999) using similar assumptions. Our model 
makes also a testable prediction: the abundance scatter 
must be smaller in the inner disk than in the outer re- 
gions, if the observed dispersion between PN abundances 
at a given galactocentric distance is due to their intrinsic 
age differences. 

In summary, we have shown that current massive star 
yields, combined with the BP99 chemical evolution model, 
reproduce fairly well most of the observed abundance gra- 
dients in the Milky Way disk. The assumptions of the 
model affect mostly points (i) and (iii) above, while the 
adopted yields affect point (ii). Some problems remain 
with C and N, and (to a smaller extent) Al. However, the 
major question of such studies, namely the evolution of 
abundance gradients, is far from being settled yet; a much 
more precise determination of ages, distances and abun- 
dances of planetary nebulae will be required for that. 
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TABLE 1. Observed and calculated abundance gradients in the Galactic disk (dlog(X/H)/dR, in dex/kpc) 



Object 
class 


Rg 
(kpc) 


He 


c 


N 


O 


Ne 


Mg 


Al 


Si 


S 


Ar 


Ref. 


HII 


4-14 


-0.001 
±0.008 




-0.090 
±0.015 


-0.070 
±0.015 










-0.010 
±0.020 


-0.060 
±0.015 


(1) 


HII 


0-10 






-0.100 
±0.020 




-0.080 
±0.020 








-0.070 
±0.020 




(2) 


HII 


0-17 






-0.111 
±0.012 












-0.079 
±0.009 




(3) 


HII 


0-12 






-0.072 
±0.006 


-0.064 
±0.009 










-0.063 
±0.006 




(4) 


HII 


5-15 








-0.040 
±0.005 














(5) 


HII 


6-9 


-0.004 
±0.005 


-0.133 
±0.022 


-0.048 
±0.017 


-0.049 
±0.017 


-0.045 
±0.017 








-0.055 
±0.017 


-0.044 
±0.030 


(6) 


B stars 


6-18 








-0.070 
±0.010 














(7) 


B stars 


5-14 




-0.035 
±0.014 


-0.078 
±0.023 


-0.067 
±0.024 




-0.082 
±0.026 


-0.045 
±0.023 


-0.107 
±0.028 






(8) 


PNI 


4-14 








-0.030 
±0.007 


-0.004 
±0.008 








-0.075 
±0.008 


-0.060 
±0.008 


(9) 


PNII 


4-14 








-0.058 
±0.007 


-0.036 
±0.010 








-0.077 
±0.011 


-0.051 
±0.010 


(10) 


PNIII 


4-14 








-0.058 
±0.008 


-0.041 
±0.008 








-0.063 
±0.010 


-0.034 
±0.010 


(9) 


All 






-0.07 
±0.02 


-0.08 
±0.02 


-0.07 
±0.01 




-0.07 
±0.01 


-0.05 
±0.02 


-0.06 
±0.01 


-0.07 
±0.02 




(11) 


Model" 
Model 6 


4-14 
4-14 


-0.003 
-0.004 


-0.062 
-0.086 


-0.093 
-0.096 


-0.061 
-0.058 


-0.082 


-0.083 


-0.090 


-0.058 


-0.053 


-0.048 





References: (1) Shaver et al. 1983; (2) Simpson et al. 1995; (3) Rudolph et al. 1996; (4) Afflerbach et al. 1997; (5) Deharveng 
et al. 2000; 

(6) Esteban et al. 1999; (7) Smartt & Rolleston 1997; (8) Gummersbach et al. 1998; (9) Maciel & Koppen 1994; (10) Maciel & 
Quireza 1999; (11) Smarrt (2000) 

Results: a Woosley and Weaver (1995) yields; b Maeder (1992) yields for He, C, N, O 



